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A  QUADRATIC-FORM  ANALYSIS  OF  THE 
COLLISIONLESS  TEARING  MODE 


I .  INTRODUCTION 

The  properties  of  the  collisionless  tearing  instability  in  the  neutral 
sheet  geometry  have  been  studied  extensively. A  physical  system  in 
which  the  instability  may  be  particularly  relevant  is  the  earth's 

magnetotall.^»^°’^^’^^  In  its  classic  form  the  inertia  of  the  current 

carriers  leads  to  the  instability.  The  equations  to  describe  the  linear 
behavior  of  the  instability  can  be  derived  using  the  standard  method  of 
characteristics.  This  yields  integro-differential  equations  for  the 
perturbed  quantities.  In  principle,  the  solution  to  these  equations  gives 
a  complete  description  of  the  linear  collisionless  tearing  instability.  In 
practice,  however,  the  equations  are  difficult  to  solve  because  evaluation 
of  complicated  orbit  integrals  is  required.  Thus,  the  majority  of  previous 
theoretical  works  have  attempted  to  extract  useful  information  without 
evaluating  the  orbit  integrals  exactly.  These  methods  include  energy 
principles,  the  use  of  approximate  orbits,  3,5,7,9,13,14 

expansion  of  orbit  integrals  in  velocity  moments. 

However,  as  the  physical  system  of  interest  becomes  more  complex,  the 

approximate  methods  of  evaluating  the  orbit  integrals  may  become  untenable. 

For  example,  if  ions  have  an  anisotropic  bi-Naxvellian  distribution  with  an 

isotropic  electron  distribution,  then  the  large  ion  orbit  contribution  to 

13 

the  instability  is  significant.  In  this  case,  the  neutral  sheet  has  to 
be  treated  by  a  3-region  (as  opposed  to  the  conventional  2-region)  matching 
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technique.  The  straight-line  orbit  and  constant-<J>  approximations  limit  the 

validity  to  relatively  small  ion  anisotropy  Tj^/T||  <  1.5.  If  the 

distribution  function  is  highly  non-Maxwell ian,  then  these  approximations 

completely  break  down  and  the  full  integro-dif ferential  equations  must  be 

solved.  In  the  problem  treated  by  Holdren^^,  the  equilibrium  orbits  and 

orbit  integrals  for  a  relativistic  neutral  sheet  were  numerically 

calculated  using  discretized  variables.  The  integro-dif ferential  equation 

for  the  perturbed  vector  potential  was  then  solved  using  an  iterative 

scheme  and  the  required  numerics  were  substantial.  More  recently,  Chen  and 

Lee^^  provided  an  integro-dif ferential  treatment  of  collisionless  tearing 

instability  in  a  highly  non-Maxwellian  neutral  sheet.  In  this  work,  the 

exact  unperturbed  particle  orbits  were  used  analytically  in  the  orbit 

18 

integrals.  Using  the  Galerkin  method,  the  integro-dif ferential  equation 

for  the  perturbed  vector  potential  was  converted  into  a  finite  dimensional 

matrix  and  solved  to  obtain  the  dispersion  relation  and  the  eigenmode 

structure.  It  was  found  that  the  eigenmode  is  strongly  localized  to  a 

1/2 

region  of  the  order  of  (P^^p)  at  the  null  plane  and  that  the  eigenmode 
is  highly  structured  in  this  region,  where  is  the  electron  Larmor  radius 
in  the  asymptotic  magnetic  field  and  x^  is  the  neutral  sheet  half¬ 
thickness.  For  the  dispersion  relation  was  found  to  have  the  general 

form  Y/kVgll  =  (p^/Xp)D  where  D  is  a  function  of  temperatures  only  and  v^^ 
is  the  electron  thermal  velocity  along  B  ,  the  equilibrium  magnetic  field. 
For  this  problem,  the  usual  straight-line  orbit  and  constant-t|« 
approximations  would  indeed  have  been  invalid.  Another  point  of  interest 
to  note  is  that  the  contribution  to  the  perturbed  current  density  J^^(x)  at 
a  given  point  x  primarily  comes  from  those  axis-crossing  orbits  with 
turning  points  at  x.  It  is  easier  to  calculate  such  contribution 
accurately  using  analytic  orbits  than  using  discretized  variables. 


In  Ref.  15,  due  to  the  relative  simplicity  of  the  geometry,  the  orbit 

integrals  and  matrix  elements  could  be  calculated  independently  of  k  and  u. 

Thus,  only  minimal  numerical  calculations  were  required.  However,  in  more 

complex  systems,  the  required  computation  may  be  more  substantial,  i.e., 

the  matrix  elements  may  have  to  be  evaluated  for  each  iteration  for 

different  values  of  u/kv  to  determine  the  dispersion  relation.  In  such 

cases,  it  would  be  desirable  to  determine  some  overall  stability  properties 

without  extensive  computation.  In  addition,  if  unstable,  it  would  be 

desirable  to  determine  the  lower  and  upper  bounds  for  the  growth  rate. 

Such  considerations  have  led  to  formulation  of  various  energy  principle 

19 

methods.  For  example,  for  ideal  MHD,  the  energy  principle  method  seeks 

to  determine  the  stability  of  a  magnetized  plasma  by  considering  the 

potential  energy  associated  with  perturbations  of  an  equilibrium.  For  the 

collisionless  tearing  mode,  a  number  of  energy  principle  analyses  have  been 
3  6  8 

performed.  ’  ’  These  methods  generally  provide  the  sufficient  condition 
for  stability.  In  these  works,  approximate  orbit  integrals  are  used,  or  in 

g 

certain  simple  cases,  no  orbit  integral  is  needed.  If,  however,  no  valid 
approximations  are  available  for  the  necessary  orbit  integrals,  then  these 
methods  do  not  yield  accurate  results  or  estimates  of  growth  rates. 

In  this  paper,  we  consider  another  method  which  is  analogous  to  energy 
principle  methods  but  which  is  more  versatile.  In  the  present  method,  we 
evaluate  the  orbit  integrals  using  the  exact  equilibrium  orbits  and  no 
trial  functions  are  used.  It  generally  yields  a  sufficient  condition  with 
a  lower  and  upper  bounds  for  the  growth  rate  for  unstable  systems. 
However,  the  method  can  also  be  iterated  to  give  accurate  growth  rates  for 
the  unstable  modes  and  the  eigenfunctions. 

The  method  depends  on  certain  general  features  of  the  eigenmodes.  One 
important  feature  is  that  the  eigenmodes  are  generally  well  localized.  The 
technique  is  potentially  applicable  to  a  wide  variety  of  systems.  The 
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basic  ideas  are  simple  and  ve  first  describe  the  method.  The  objective  of 
this  paper  is  to  provide  a  rigorous  theoretical  basis  for  the  method.  For 
this  purpose,  ve  prove  a  number  of  general  properties  of  the  quadratic  form 
based  on  a  self-adjoint  integro-dif ferential  operator.  Ve  then  illustrate 
the  technique  by  applying  it  to  the  collisionless  tearing  mode  in  a  highly 
non-Maxvellian  neutral  sheet  that  has  been  solved  previously. 

II.  A  Quadratic  Form  Formulation 

In  this  section,  ve  present  the  formalism  for  a  quadratic-form  method 
for  treating  the  collisionless  tearing  instabilities.  Ve  first  describe 
the  general  framevork  and  then  provide  a  specific  application.  For  the 
sake  of  concreteness,  ve  vill  use  a  neutral  sheet  geometry.  The  method  can 
be  applied  to  other  configurations. 

Figure  1  shovs  the  magnetic  field  and  the  coordinate  system.  The 
equilibrium  magnetic  field  is  given  by  B  (x)  =  B  (x)  z  vith  B  (-x)  =  -B(x) 

O  2  “  2 

* 

and  B,  =  0  at  X  =  0.  The  current  density  is  J_(x)  =  J  (x)y  vith  J  (-x)  = 
J^(x).  Because  of  pressure  balance,  the  particle  density  n^(x)  and  current 
density  Jq(x)  are  peaked  about  x  =  0  vith  a  characteristic  half-thickness 
of  Xp.  For  |x  I  >  Xp,  B^  is  nearly  constant  and  the  current  density  is  much 
smaller  than  for  |x|  <  Xp.  The  tearing  instability  in  a  neutral  sheet  can 
be  described  by  perturbations  of  the  form  <|>(x,z,t)  =  tlrCxjexpCikz  -  iwt), 
vhere  the  vave  vector  k  =  k  z  is  taken  to  be  parallel  to  the  equilibrium 
magnetic  field.  Ve  assume  that  all  perturbations  vanish  as  t  -<*>.  In 
this  paper,  ve  vill  neglect  the  scalar  potential.  The  inclusion  of  the 
scalar  potential  gives  rise  to  another  coupled  integro-dif ferential 
equation  vhose  basic  character  is  similar  to  that  of  the  equation  to  be 
treated  here.  Thus,  the  present  method  can  be  generalized  to  include  the 
scalar  potential  in  a  straightforvard  manner.  In  reality,  the  neglect  of 
the  scalar  potential  implies  that  high  perturbation  frequencies  (|w|  > 


is  the  ion  cyclotron  frequency  in 


cannot  be  treated  adequately,  where 
the  asymptotic  magnetic  field. 

The  Ampere's  lav  gives 


^  -  khw  +  .  0  ,  (1) 

dx  ^ 

where  the  perturbed  vector  potential  is  ■  ^(x,z,t)}r  and  ^ly  is  the 

perturbed  current  density.  It  is  in  evaluating  that  one  must  evaluate 
orbit-integrals.  These  integrals  may  be  regarded  as  linear  operators 
acting  on  4).  Thus  we  can  re-write  Eq.  (1)  as 


L«|*  =  0, 


(2) 


where  L  is  an  integro-differential  operator  to  be  determined  for  each 
specific  problem.  We  can  also  define  a  quadratic  form  I  by 

K  +  IO/)  ■  I  dx^*(x)L^(x)  ,  (3) 

where  ^  is  the  complex  conjugate  of  It  is  clear  that  Eq.  (2)  can  be 

recovered  by  the  variational  principle  requirement  »  0,  which 

automatically  implies  I('|'|*)  =  0.  If  we  let  be  an  eigenfunction  of  L 
corresponding  to  the  eigenvalue  then  we  can  write 


where  C  is  the  linear  operator  describing  the  orbit  integrals.  Here,  we 
have  used  0  as  |x|  ->  «>  and  assumed  ^  to  be  normalized  to  unity.  Note 

that  the  eignevalues  depend  on  k  and  a  and  that  they  are  primarily 
determined  by  the  first  and  third  terms  of  the  above  equation.  However, 
Eq.  (2)  can  also  be  thought  of  as  an  eigenvalue  equation  L<|/(x)  =  Xg/(x), 

where  X  =  X(k,  (o).  Equation  (2)  is  then  equivalent  to  setting  X  -  0. 
Thus,  the  dispersion  relation  can  be  obtained  from 

X(k,  «)  =  0.  (5) 

In  this  paper,  we  will  prove  a  number  of  mathematical  properties  of  the 
eigenvalues  X^  and  the  integro-differential  equation.  These  properties 
will  be  used  in  a  method  to  determine  the  stability  condition  for  a  given 
system  with  respect  to  the  tearing  mode.  It  is  an  objective  of  this  paper 
to  provide  a  sound  theoretical  basis  for  the  method  and  to  illustrate  the 
application. 

Following  the  treatment  of  Ref.  15,  we  use  an  equilibrium  distribution 


function  of  the  form  Fj^  =  ^^llj ^^||j ^ '  where  j  =  e,  i  and 

(Hj^j,Pyj,  H||j)  are  the  single-particle  constants  of  motion.  Here,  Hj^j  = 

2  2  0  0 

(mj/2)(v^  +  Vy  ),  Pyj  =  m^Vy  +  (qj/c)Ay(x)  where  Ay(x)  is  the  equilibrium 


vector  and  H||j  =  (mj/2) 

Maxwellian  with 


For  simplicity,  we  will  choose  F,,  to  be 


F„j(H„)  =  (2nT„/mj)"^^^exp(-H„/T„).  (6) 

For  the  species  j,  m^  is  the  mass,  qj  is  the  charge  and  Vj  is  the  mean 
drift  velocity  in  the  y  direction.  In  the  remainder  of  the  paper,  we  will 
omit  the  species  index  j  when  no  confusion  results.  Then,  it  can  be 


shown'*^^  that  the  left-hand  side  of  Eq.  (1)  can  be  written  as  an  integro- 


differential  operator  L  given  by 


Ane^v  1 


E”^  [dx' [^Kj^(x,x  )  +  K2(x,x  )]  \Kx' )  . 
.  m.  c  '  111 
J  j 


Here,  the  kernels  are  given  by 


'  V  (x)v  (X  )  j. 

K,(x,x  )  =  E  t/H  'p  ^ ^ — > —  exp  -inQ(t-t  ) 

1  n  Py>  n  1  L  J 


>  =  ^Jt(h,,p')(^  ^  i^2<‘’n>  357  f  7,  'r\'i  >  <9> 

n  '  1’  y'  jll  i  |v^(x)v^  (x  )|  ■- 


where  the  orbital  period  is  given  by 


■  f  TWT  ’ 


*  («  +  nQ)/kVjll 


Here,  the  integration  is  carried  out  over  one  complete  cycle.  The 


coordinates  (x  ,  v  )  represent  the  particle  orbits  in  the  equilibrium  field 


with  the  conditions  x  (t  =  t)  =  x  and  v  (t  =  t)  =  v.  The  dispersion 


functions  are  given  by 


Z(0  =  J  dt  ^ 


(A 

j^f 


E 


ir 

’v»’ 


E 


B5 


Mil)  =  1  +  CZ(t). 


In  the  above  expressions,  is  the  thermal  velocity  of  the  species  j 
associated  with  the  motion  parallel  to  the  equilibrium  magnetic  field.  In 
deriving  Eqs.  (7),  (8)  and  (9),  we  have  used 


Vjj(x  )»J»(x  )  =  ^  ♦^(x)exp  {inQ(t  -  t  )1 


f  I  )Iexp[-inQ(t  -  x)] 


where  T  =  T(Hj^,  P  )  and 


Q  =  2ii/T(H^,  Py), 


with  the  orbital  period  T  given  by  Eq.  (10). 


In  writing  Eq.  (11),  use  has  been  made  of  the  fact  that  an  equilibrium 

r 

orbit  can  be  thought  of  as  a  function  of  the  present  coordinates  at  t  =  t. 
The  kernels  and  K2  generally  must  be  evaluated  numerically.  However, 
there  is  no  need  to  follow  equilibrium  orbits.  The  velocity  components  can 
be  evaluated  analytically  as  functions  of  H.  and  P  .  In  doing  the  velocity 

y 

space  integration,  all  the  equilibrium  orbits  are  included  exactly. 

It  is  easy  to  see  that  the  third  term  on  the  right  hand  side  of  Eq.  (7) 
is  the  so-called  adiabatic  term.  The  third  and  fourth  terms  combined  give 
the  perturbed  current  density.  In  deriving  Eqs.  (8)  and  (9),  we  have 


•'1. 


included  all  harmonics  of  the  orbital  periods.  However,  the  conventional 


tearing  mode  corresponds  to  the  n  =  0  case.  In  this  paper,  we  will  also 
set  n  =»  0  henceforth  for  simplicity  and  restrict  ourselves  to  the  purely 
growing  mode,  »  =  iy. 


III.  A  Non-Maxwellian  Neutral  Sheet 

If  we  set  n  =  0  and  w  =  iy.  the  operator  L  is  self-adjoint  since  it  is 
real  and  symmetric  in  x  and  x  .  Then,  its  eigenvalues  are  real  and 
eigenfunctions  can  be  chosen  to  be  real  and  orthonormal  without  loss  of 
generality.  Let  ♦jjCx)  be  such  an  eigenfunction  with  the  eigenvalue  X^. 
Then, 


m  m  m 


k2  ^  ^ 

c 


i;s.Jdx*2(x)jd\vy  jjji 


+  Jdxdx  ^  K2(x,x  ),  (13) 

j  m^  c  llj 

where  we  have  used  0  as  |x|  -»  ®.  Note  that  the  dependence  on  m 

occurs  only  through  the  eigenfunctions 

We  will  demonstrate  the  following  properties.  We  will  first  show  that 

the  eigenvalue  X^  increases  without  bounds  as  y  approaches  for  all  m  and 

k.  Then,  we  will  show  that  X^  approaches  an  asumptotic  value  for  each  k 

and  m  as  y  approaches  +®  and  that  the  asymptotic  values  are  bounded  from 

above  but  not  from  below.  We  then  prove  that,  for  any  finite  value  of  y, 

there  exist  only  a  finite  number  of  m  values  for  which  X  is  positive. 

m 


.1  •  .1  ■  .>.  .1.  ,t. 

::: 


We  now  prove  that  X  ■*+<*>  for  all  m  as  y  First,  we  consider  an 

Dl 


analytic  continuation  of  Eq.  (13)  to  y  <  0.  The  first  three  terms  of  Eq. 


(13)  have  no  explicit  dependence  on  y.  Equation  (9)  shows  that  the  <2 
also  has  no  y  dependence  if  n  =  0.  Equation  (8)  shows  that  depends  on  y 


through  the  dispersion  function  V.  For  large  negative  y,  we  have  the 


asymptotic  expression 


w(^)~  -  2ii^^^y  exp(Y^). 


Therefore,  as  y  we  have  W  ^  In  addition,  can  be  re-written 


<^(x,x  )  =  W(b^)  j:i^[jdx+^(x) 


where  b^  =  «/kVj||.  Since  T  is  positive  definite,  the  integral  is  positive 

t 

definite.  Thus,  we  see  Kj^(x,x  )  ->  +  ®  in  the  limit  y  -  ®.  Because  all 


other  terms  in  Eq.  (13)  are  independent  of  y,  we  conclude  that 


lim  X  =  +® 
™ 

Y->-«D 


for  all  m.  In  the  other  limit  y  +  ®,  we  see  that  Z(b^)  -*  0  and  W(bp)  -» 


-i-l  so  that,  for  a  given  k,  X^  reaches  an  asymptotic  value  for  each  m.  In 


practice,  the  dependence  of  y/kv  becomes  unimportant  for  y/kv  >  4  where 


Z(4i)  =  0.1.  The  large-Y  asymptotic  expression  for  X^  then  becomes 
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It  is  clear  that  the  asymptotic  values  of  may  be  positive  or  negative. 

It  is  of  interest  to  note  that  if  and  3Fj^/dHj^  are  integrable  at  least  in 

the  ^-function  sense,  then  the  asymptotic  values  of  X  are  bounded  from 

n 

above.  However,  they  need  not  be  bounded  from  below  since  it  is  possible 

2 

to  make  the  integral  of  (d^^/dx)  arbitrarily  large  by  introducing 
"wiggles".  Since  lim  X_  ■*  *<»  for  all  m,  and  X_  can  be  negative  for  some  m 

^  in  JD 

as  r  -»  -f  a,  it  is  clear  that  some  eigenvalues  must  cross  zero.  If  the 

zero-crossings  occur  for  y  >  0,  then  the  eigenfunctions  corresponding  to  X 

«  0  provide  the  solutions  to  the  tearing  mode  equation  (2).  The  values  of 

Y  where  the  crossings  occur  are  the  corresponding  growth  rates. 

We  now  demonstrate  an  important  property  that  for  any  finite  y, 

negative  or  positive,  there  are  only  finitely  many  positive  eigenvalues  X^^. 

For  any  given  k  and  a  finite  Yi  it  is  easy  to  see  that  only  the  first  term 

of  Eq.  (13)  can  be  infinite.  Thus,  it  is  sufficient  to  show  that  there  are 

2 

only  finitely  many  m's  for  which  the  quantity  Jdx(d*|^/dx)  is  finite. 

Let  {U|^(x))  be  a  complete  set  of  basis  functions.  Then,  we  can  write 


k>l 

It  is  Important  to  consider  the  nature  of  the  eigenfunctions.  We  note  that 

2 

Eq.  (1)  can  be  cast  in  the  form  of  a  Schroedinger  equation  with  -k  playing 
the  role  of  the  energy  of  a  particle  in  a  potential  well  -(Aii/c)!^^.  For 


the  tearing  mode  problem  at  hand,  the  perturbed  current  is  sharply  peaked 

at  the  magnetic  null  plane  (x  >  0).  Thus,  ve  seek  bound  state  solutions  of 

Bq.  (1)  localized  around  x  »  0.  As  a  general  remark,  ve  note  that  the 

equilibrium  distributions  are  such  that  the  current  and  particle  densities 

are  localized  near  x  0  with  a  characteristic  half-thickness  x  as 

P 

discussed  earlier.  Then,  in  the  operator  L,  only  the  first  two  terms  are 
important  [Eq.  (7)].  For  |x|  »  x^,  the  current  density  is  essentially 

zero  and  the  perturbed  eigenfunction  vanishes  exponentially.  As  a  simple 
model  with  these  features,  ve  use^^ 

.  (iiuiQ/2ll)  (H  -  VjPy  -  Tj^) 

vhere  and  Tj  are  constants  and  nQ  »  0^(0).  This  distribution  is  such 
that  nQ(x)  and  Jq(x)  are  constant  for  |x |  <  x^  and  zero  for  |x|  >  x^  for 

each  species,  vhere 

P  ”  4llen^(Vj  -  Vg)(qVj)/c2 


Here,  q  is  the  charge  of  the  species.  For  this  distribution  function,  Eq. 
(1)  can  be  solved  exactly  outside  the  plasma  sheet  (|x|  >  x^)  to  give 


♦  =  ♦(Xp)expl-k( |x|-Xp)]. 


Thus,  ve  can  vrite 


P»  .2  -’'p  .2  , 

J  '•’■(sr)  ■  I  '"‘t)  * 


-X 


Ve  define 


•  I 


(17) 


-X 
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Then,  is  finite  if  and  only  if  is  finite.  For  simplicity,  we  treat 

♦  as  if  it  were  entirely  confined  in  the  region  (-x  ,  x  ).  Then,  the 
m  P  P 

Fourier  components  for  even  are  simply 


Uj^(x)  .  d  ^^^cos^^j. 


(18) 


Substituting  Eq.  (18)  into  Eq.  (17),  we  obtain 


m  (d  j  ^  “mk 
k-1 


(19) 


The  orthonormality  of  the  eigenfunctions  4  requires 

B 


E  V  “ 


(20) 


k.l 


for  all  a.  In  addition,  we  have 

E  ■  *•  (21) 

m^l 

for  all  k.  This  can  be  seen  by  noting  that  the  matrix  A  ■  is  such 

T 

that  AA  -  I  from  the  orthonormality  of  eigenfunctions  Since  A  is  a 

T 

square  matrix,  A  is  the  inverse  of  A,  i.e.,  A  is  an  orthogonal  matrix. 

•} 

Then,  Eq.  (21)  follows.  From  Eqs.  (19)  and  (20),  we  see  >  (n/d)  for 

2 

all  m  with  R^^  -  (R/d)  if  and  only  if  a^i^j  =  1  and  a^^j^  •  0  for  all  k  ^  1. 

We  now  prove  that  there  exist  only  finitely  many  m  for  which  R^^  is 

Define  the  following 


finite.  Consider  a  set  S  of  all  finite  R  . 

HI 

quantities 


and 


"  k=N+l 


N 


so  that  +  Q^.  Because  the  sequence  R^  converges  to  R_^  as  N  -»  », 


we  see  that  for  any  €  s  and  any  arbitrary  e  >  0,  there  exists  a  finite 
integer  N  such  that  Q„  <  c  for  all  N  >  N  .  Now,  consider  a  finite  subset 
of  S  such  that  has  M  elements.  Let  be  the  average  value  of  the 

elements  of  Sj^.  Then,  by  denoting  the  elements  of  by  m  =  1,...  M,  we 
have 


<Vm 


M 

m=l 


H  N  M 

=  M  ^  M  ^  ^N* 

m=l  k=l  m=l 

Here,  for  any  e  >  0,  N  is  chosen  such  that  Q„  is  less  than  e  for  all  R  in 

w  m 

S^.  Interchanging  m-  and  k-  summations  and  summing  over  m  to  infinity,  we 
obtain 

N  ® 

<R„>H  <  R  E  E  .  <0„>„ 

k=l  m=l 

where  is  the  average  over  Then,  using  Eq.  (21)  and  carrying  out 

the  k-summation,  we  have 


«„>«  <  ^  <Vh 


(22) 
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for  all  Rjji  in  Sj^.  Clearly  for  any  c,  it  is  easy  to  show  <Qjj>^  <  c  by 

2 

induction.  We  can  choose  e  <  (n/d)  .  Nov,  we  suppose  S  has  infinitely 

many  elements.  Then,  we  can  take  the  limit  M  -»  ».  Recall  that  R  >  (n/d)^ 

m  ~ 

for  all  m. 

2 

Therefore,  lim  <R_>m  >  (n/d)  .  However,  the  inequality  (22)  shows  that 

lim  <  c  <  (n/d)^ 

M  -»  ® 

This  is  a  contradiction  due  to  the  supposition  that  N  is  infinite. 

Therefore,  S  cannot  have  infinitely  many  elements. 

The  properties  of  described  above  suggest  a  practical  method  for 

assessing  stability  of  a  neutral  sheet  when  detailed  computation  is  time- 

consuming.  First,  compute  the  positive  eigenvalue  spectrum  at  y  -  0. 

Suppose  there  are  n^  positive  eigenvalues.  Next,  compute  the  positive 

eigenvalue  spectrum  at  a  large  y  «  Y2  with  n2  positive  eignevalues.  If  n^ 

and  n2  are  not  equal,  then  there  are  at  least  |n2  -  n^|  values  of  y  >  0  for 

which  Eq.  (2)  is  satisfied.  That  is,  there  are  at  least  |n2  -  n^ |  unstable 

modes.  This  method  would  be  particularly  useful  when  evaluation  of  kernels 

for  each  value  of  k  and  (o  is  numerically  prohibitive.  In  cases  where  some 

numerical  computation  is  still  practical,  the  above  procedure  can  be 

Iterated  by  considering  successively  smaller  intervals  Ay  ■  (y2  -  y^)* 

From  Eqs.  (3)  and  (17),  we  see  that  there  can  be  only  a  finite  number  of  m 

such  that  is  finite  for  finite  y.  In  particular,  there  can  be  only 

finitely  many  positive  X^^  for  y  =  0.  There  are  infinitely  many  m  such  that 

X  's  cross  zero  for  y  <  0.  That  is,  there  are  infinitely  many  stable  modes 
in 

and  only  a  finite  number  of  stable  modes.  It  is  also  clear  that  if  the 
dimensionality  of  the  matrix  is  large  enough,  the  values  of  the  positive 
eigenvalues  are  not  significantly  affected  by  the  dimensionality  of  the 
matrix.  That  is,  we  only  need  to  solve  a  finite  dimensional  eigenvalue 


equation.  Thus,  Eq.  (2)  can  be  reduced  to  a  finite  dimensional  matrix 
equation  to  find  all  the  unstable  modes,  if  the  system  is  unstable.  If  Ay 
can  be  made  sufficiently  small,  then  we  obtain  the  growth  rate  and  the 
eigenmode  structure  accurately. 

IV.  Summary  and  Discussion 

We  have  described  a  quadratic-form  method  which  may  be  useful  for 
studying  the  linear  collisionless  tearing  mode  stability  of  a  neutral  sheet 
in  cases  where  the  complicated  orbit  integrals  make  it  impractical  to 
calculate  the  eigenmode  structures.  Ve  have  provided  the  necessary 
mathematical  justification.  This  method  may  be  Iterated,  where  practical, 
to  provide  a  necessary  and  sufficient  condition  regarding  the  stability  of 
a  complicated  system  and,  if  unstable,  the  growth  rates. 

In  this  paper,  we  have  used  a  mathematically  tractable  distribution 

function.  However,  the  basic  conclusions  of  the  paper  are  not  limited  to 

the  Idealized  model.  In  the  example  given  above,  the  kernels  and  K2 

depend  on  the  standard  dispersion  functions  Z  and  V  because  of  the 

Haxwellian  distribution  in  H||.  As  a  result,  the  asymptotic  behavior  of 

and  K2  for  y  ±  *  can  be  determined  readily.  Similar  asymptotic  behavior 

occurs  in  other  systems  in  which  the  distribution  functions  are  not 

Maxwellian  in  H||.  It  has  been  shown  that  the  behavior  of  the  eigenvalues 

X  is  determined  by  the  quantity  R  [Eq.  (17)].  Ve  have  shown  that  R  is 
n  In  n 

finite  only  for  a  finite  number  of  m.  This  behavior  is  due  to  the  fact 
that  eigenfunctions  A  are  localizable  to  a  finite  region,  which  in  this 
case,  is  the  magnetic  null  region.  Thus,  the  properties  of  which 

are  necessary  for  the  method  can  be  generalized  to  a  wide  variety  of 
systems.  Ve  expect  that  the  present  quadratic-form  technique  can  be  used 
to  analyze  more  realistic  physical  systems.  In  this  paper,  ve  have  shown 
how  the  necessary  properties  must  be  assessed  for  application  to  a  given 


system. 
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